Generation of two-color polarization-entangled optical beams with a self-phase-locked 

two-crystal Optical Parametric Oscillator 
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A new device to generate polarization-entangled light in the continuous variable regime is intro- 
duced. It consists of an Optical Parametric Oscillator with two type-II phase-matched non-linear 
crystals orthogonally oriented, associated with birefringent elements for adjustable linear coupling. 
We give in this paper a theoretical study of its classical and quantum properties. It is shown that 
two optical beams with adjustable frequencies and well-defined polarization can be emitted. The 
Stokes parameters of the two beams are entangled. The principal advantage of this setup is the 
possibility to directly generate polarization entangled light without the need of mixing four modes 
on beam splitters as required in current experimental setups. This device opens new directions 
for the study of light-matter interfaces and generation of multimode non-classical light and higher 
dimensional phase space. 

PACS numbers: 03.67.Hk, 03.67.Mn, 42.65.Yj, 42.50.Dv 



I. INTRODUCTION 

Developing quantum memories constitutes an essential 
milestone on the route towards quantum communication 
networks. As it is difficult to directly store photons, the 
quantum information has to be stored in a material-based 
quantum system. A particularly promising application 
for polarization-entangled light is the ability to couple it 
with atomic ensembles: the algebra describing the quan- 
tum properties of polarized light via Stokes operators 
ELH>U is exactly the same as that describing the quan- 
tum properties of atomic spin be it single or collective 
spins. Quantum state exchange between light fields and 
matter systems has been recently experimentally demon- 
strated 0, Q and it has been shown that such systems 
can be used as quantum memories |(|. In addition, from 
a technical point of view, the use of polarization states 
offers simpler detection schemes, without the need for 
local oscillators. These features make the study and ex- 
perimental generation of non-classical polarization states 
of first importance for continuous variable quantum com- 
munication. 

Furthermore, the study of higher dimensional phase 
space is a fascinating and promising subject. Indeed, it 
has been shown that multimode non-classical light could 
be a very powerful tool for an efficient processing and 
transport of quantum information Q. Already exper- 
imental implementations and applications of such non- 
classical multimode light have been demonstrated in the 
continuous- wave regime 0, 0, 0, [n| . Such systems ci- 
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ther involve a very large number of modes [ToL lllj or a 
much smaller number (three in the case of ref. |9(). In- 
creasing the number of correlated modes is of particular 
importance, for instance in the case of optical read-out 
technologies ^2 an( i super-resolution techniques 01 ■ In 
this context, generating continuous- variable polarization- 
entangled light ' which exhibits three or four mode 
correlations, is a first step towards increasing the phase 
space dimension. 

Such polarization-entangled beams in the continuous 
variable regime have been recently produced experimen- 
tally with injected type-I phase-matched Optical Para- 
metric Oscillators (OPOs) below threshold [pi Hfjj or 
using x^ 3 ) effects in a cloud of cold atoms 0. All these 
methods require linear interferences of two quadrature- 
entangled modes with two bright coherent states, which 
limits the amount of entanglement reached and the sim- 
plicity and scalability of the setup. 

We investigate here a new approach to directly gener- 
ate polarization-entangled light without the need of mode 
mixing. This approach stems from an original device we 
previously studied, a self-phase-locked optical paramet- 
ric oscillator. This system consists of a type-II phase- 
matched nonlinear x*- 2 -* crystal in an optical cavity and 
emits fields that are phase-locked through linear coupling 
[PH ITcl H^: a birefringent plate, which can be rotated 
relative to the principal axis of the crystal, adds this 
coupling between the orthogonally polarized signal and 
idler beams and results in a phase-locking phenomenon 
that is well-known for coupled mechanical or electrical 
oscillators [2(J. Such systems have recently attracted a 
lot of attention as efficient sources of non-classical light 
pH [52, ■ Stable operation has been demonstrated ex- 
perimentally even with very small coupling [ItI IziL l25j 
and their non-classical properties are very encouraging 
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FIG. 1: Left: Ring cavity two-crystal OPO : c a ,p are type-II nonlinear crystals, are A/2 waveplates at ^ and A waveplates 
at loq. Mi,3 are highly reflective mirrors at u> ai t, and M2 is a partially reflective mirror at uj a ,b- All mirrors are transparent at 
ojq. G is a diffraction grating which separates the two frequency components. Each of them have two orthogonally-polarized 
bright components. Center: neutral axes seen by the fundamental wave at frequency cjq. Right: neutral axes seen by the 
sub-harmonic waves at frequencies uj\ and u>2. Dotted lines denote fast axes and continuous lines slow axes. 



|24| . We propose here a new implementation where two 
crystals, and not only one, arc placed with their neutral 
axis orthogonal in a cavity containing two bircfringent 
waveplates for adjustable linear coupling. Each crystal 
generates a pair of signal and idler beams and the wave- 
plates couple together the two signal beams and the two 
idler beams so that phase locking can be achieved in a 
non-frequency degenerate operation. As we will show be- 
low, this device allows for a direct and efficient generation 
of polarization-entangled optical beams. 

In this paper, we investigate in detail the classical and 
quantum properties of the system. The paper is orga- 
nized as follows. In Sec. ^ we begin by presenting the 
linear and non- linear elements of the two-crystal OPO. A 
propagating Jones matrix is associated to each elements 
and permits to access the round-trip matrix and the 
stationary solutions. Section ITTT1 considers the quantum 
properties of the emitted beams in terms of quadrature 
operators in the Fourier domain. The usual input-output 
linearization technique is used and quantum correlation 
spectra are formulated. In Section llVI these quadrature 
correlations are then interpreted in terms of polarization- 
entanglement through the Stokes operators. Considering 
various criteria, we demonstrate theoretically that this 
device is an efficient source of polarization-entanglement. 

II. SETUP AND CLASSICAL PROPERTIES 

In this section, we present the basic scheme of the two- 
crystal OPO. Round-trip equations and stationary solu- 
tions are then derived by using Jones matrix formalism. 

A. Linear and nonlinear elements 
in the OPO cavity 

The setup is sketched in Fig. ^ Two identical type-II 
phase matched (x' 2 ') crystals oriented at 90° and two 
bircfringent waveplates are inserted inside a ring cavity. 
Signal and idler fields are resonant but the pump is not 
enhanced. The orientation of the various axes with re- 



spect to one another are critical. They are summarized 
on Fig. n 

• the nonlinear crystals c a and cp have their fast axes 
orthogonal; 

• waveplate l\ has its fast axis parallel to the fast axis 
of the crystal c a ; 

• the fast axis of the waveplate I2 makes an angle p 
with respect to the fast axis of crystal cp. In this 
paper, we will restrict ourselves to small values of 
P- 

Four fields may propagate inside the cavity: their clas- 
sical field amplitudes are denoted ai ; 2 and b\ y 2 where the 
index corresponds to the polarization and the letter to 
the frequency u a or cut, (Fig. The pump field is po- 
larized at 45° of the crystals' axes and its amplitude is 
denoted ao- We consider that the pump and subharmonic 
fields are coupled via the nonlinear crystals. In c a , the 

(x) 

phase matching is such that a is coupled to a\ and 62 
while in cp, is coupled to and b\. 

We suppose that the waveplates for the sub-harmonic 
waves have identical dephasing for the two fields. M\ 
and A/3 arc highly reflective mirrors while the amplitude 
reflection coefficient for M2 is denoted r. For the sake 
of simplicity, this coefficient is assumed to be real and 
independent of the frequency and of the polarization. 



B. Propagation matrices 

Bircfringent elements are usually described by Jones 
matrices |2(j. The total field can be decomposed into 
four components corresponding to the two sub-harmonic 
frequencies and the two linear polarizations. The wave- 
plates l\ and I2 couple the fields a\ and 02, and b\ and b\ 
independently while the crystals couple a\ and b^, and 
d2 and b\ independently. Thus, the basis {di, b\, (12, b\} 
is sufficient. 
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The waveplate l\ is a A/2 with its axes parallel to the 
basis axes so that its associated matrix is diagonal: 



M h = 



\ 





—ik^ne 






l\ ' ' k i I'll 

— te b 



(i) 



where k a ,b denotes the wave vectors, n the mean index 
of refraction (dispersion is neglected) and e the thickness 
of the plate. 

To determine the transfer matrix of the nonlinear crys- 
tals, the phase-matching will be taken perfect. We recall 
that the field amplitudes at the output of the first crystal 
can be written 



(out) 



(in) 



(x) 

ga 



and (&<"">)* = (bt ) y+9(a^y(a^ ) ) (2) 



(x) 

to the first order in g where a denotes the pump am- 
plitude at the output of the first crystal, 
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(3) 



and where g is the nonlinear coupling coefficient given by 
9 = lX {2) ' 



2c 2 eononiri2 



(4) 



Equations are only valid close to the oscillation thresh- 
old where the pump depletion is small, which corresponds 
to the usual operation regime of such non-linear devices 
in quantum optics |24j . The transfer matrix of the first 
crystal takes thus the following form: 
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ik a n,±l 
°>o r — ik h Ti2l 



^2 



V 



-ik b n 2 l 

o 
o 





a i k a 7i2 ' 

o 









\ 



(5) 

As stated before, perfect phase-matching has been as- 
sumed, ni 2 stand for the refractive indices along the 
neutral axes and I is the crystal thickness. In the second 
crystal, only 0,2 and 6* are coupled to the pump. Assum- 
ing that the two crystals are at the same temperature, 
the transfer matrix is given by: 



M c 





-ik^nil 









ik a n^l 

ga { y) e- lk » n21 






-ik b n 2 l 



Ay) 



(6) 

denotes the pump amplitude along y between the 



two crystals, 



Ay) 



V2 



(lo- 



ir) 



Finally, the last waveplate is a A/2 waveplate rotated 
by an angle (p + n/2) with respect to l±. To the first 
order in p , the transfer matrix is found to be 0] 



Mi. 



I —ie 
\ 



ik a ne 





ie e ik " ne 








ie -ik b ne 





-ie e- tk " ne 


ik a ne 





ie*'"™ 








-ie e-' k >> ne 





-ie ik » ne 



\ 



J 

(8) 

where eo = sin2p ss 2 p. 

As mentioned earlier, we assume that the mirror prop- 
erties are independent both of the frequency and of the 
polarization. The free propagation matrix is thus 



( * 



V 














-ikfrL 



(9) 



where L is the cavity length without taking into account 
the birefringent elements (waveplates and crystal). 
One can now calculate the round-trip matrix: 



M prop .M h .M Cfj .M Calpha 



■M, 



(10) 



Rather than giving the round-trip matrix in the general 
case, let us now make the following assumptions which 
are usually verified experimentally: 

• the finesse for the subharmonic fields is high. The 
coefficient r is close to 1 and we put r = 1 — n with 

• we assume that the double resonance condition is 
verified, that is both 



A a , 6 = k aM (L + 2ne + l(ni + n 2 )) 



(11) 



are close to an integer multiple of 27r. We denote 
5 a , b their small round-trip phase dctunings. 

• the waveplates l\ and I2 are A/2 for the two fre- 
quencies u a and Lib- 
ia this case, the round-trip matrix to the first order in 
5 a ,b,g, k and eg takes the following form: 



Mr 



/ 1 + i5 a - K g^L 

gf 2 l-ifc- 
-eo 
e 



K 

1 + i5 a — K 

9 V2 C 







an 

- iSb 



■iip 



(12) 

where ip = (kbn\ + k a ri2)l- The phase of the pump 
amplitude between the two crystals has been shifted by 
(k a + h)ne. 



C. Stationary solutions 

The stationarity of the solutions corresponds to the 
fact that the field amplitudes after one round trip must be 
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equal to the initial ones. This condition can be expressed 
in terms of the round-trip matrix: 



M r 



ai \ 








- 


h* 








h\ ) 







(13) 



This system has a non-trivial solution only if 

det(M rt - 14) = (14) 
which leads to a stringent condition on the dephasings 

5 a = 5 h = 5 (15) 
and provides also an expression for the pump power 

T (threshold) 1 1 2 

1 o = \ a o\ ■ 

We will restrict ourselves to the case where ip = 0[27r], 
which corresponds to the lowest threshold area. Eg 1141 
yields to: 

j-(threshold) 



\ a o\ 



{ K 2 + (5±e ) 2 ). (16) 



This equation has two solutions, corresponding to two 
possible regimes. In this paper, we will focus on the 
regime of lower threshold. 

These equations are similar to that obtained in a one- 
crystal self-phase locked OPO, which has been studied in 
detail in |l8l Il9l I25I ] . For a given input pump intensity, 
the phase-locked operation can be obtained only for a 
given range of the relevant parameters (i.e cavity length 
and crystal temperature) which defines a so-called " lock- 
ing zone". This locking zone is typically a cavity reso- 
nance width large in terms of length and a few Kelvins 
in terms of temperature, well within experimental capac- 
ities. For a one-crystal self-phase locked OPO, phase- 
locked operation was indeed observed even for extremely 
small angles of the wave plate Q where it does not per- 
turb the quantum properties of the system |28j | . 

Let us consider the lowest threshold given by Eq. ITSl 
and reached for 5 = cq. This working point, associated 
with ip — 0[27r], can be obtained experimentally by simul- 
taneously adjusting both the crystal temperatures and 
the frequency of the pump laser. Equation 1161 becomes 
thus: 

j-(threshold) 







(17) 



M r 



At this working point, the round-trip matrix is simpler 
and can be rewritten: 

/ 1 + ieo — k k eo \ 

K 1 — ieo — K — eo 

— eo 1 + ieo — k K 

eo K 1 — ieo ~ K / 

(18) 

The eigenvector of the round-trip equation, M rt J = 
J ', is given by 



V 



/ 1 



J = J 



1 

—i 
—i 



(19) 



where J is a (a priori complex) constant. This expres- 
sion shows that the field generated consists in two modes 
at frequencies ui a and u>b with right circularly polariza- 
tion. J is determined by the pump depletion equation: 



ao = (a™ - gaib 2 ) e 



i kilo I 



(20) 



where a™ is the input pump amplitude. In order to solve 
this equation, one can set the pump phase between the 
two crystals to be zero : this amounts to changing the 
phase origin which will not change the properties of the 
system. The previous equation yields 



= 1 



(in) 



(21) 



with 4 m) 

tion 



(in) 



One then gets the usual OPO solu- 



J = 



K 

9 



l)e 



i<t>i 



(22) 



where a reduced pumping parameter a has been defined 
equal to the input pump amplitude normalized to the 
threshold 



\ 



Jin) 



^threshold) 



(23) 



When a is larger than 1, i.e. when the OPO is pumped 
above a defined threshold, four bright beams with fixed 
relative phases given by Eq. ^|are generated. Since the 
two signal (resp. idler) beams have a fixed phase relation 
and identical frequencies, they form a single beam with 
a circular polarization state. The system thus generates 
two circularly polarized beams which arc frequency tun- 
able. 

The next sections are devoted to the study of the quan- 
tum properties. Correlations and anticorrelations of the 
couple signal/idler emitted by each crystal, i.e. (01,62) 
or (a2,£>i), are formulated. 



III. QUANTUM PROPERTIES AND 
QUADRATURE FLUCTUATIONS 

In order to calculate the fluctuations for the involved 
fields when the system is pumped above threshold, we 
appl y th e usual input-output linearization technique (see 
e.g. [27j )• Individual noise spectra as well as the corre- 
lations between any linear combinations of the various 
fluctuations are derived in the Fourier domain. 



A. Linearized equations 

We first determine the classical stationary solutions 
of the evolution equations : these solutions are denoted 
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a[ s \b^ where i = 0,1,2; j — 1,2. We then linearize The evolution equations for the fluctuations arc de- 

the evolution equations around these stationary values duced from the round-trip matrix taking into account 

by setting : a, = + 8a u i = 0, 1, 2 and bj = + th c stationary solutions [131122 
Sbj,j = 0, 1. 



dSa\ 
dt 

dSa2 
dt 

d5b 1 
dt 
_ d5b 2 
dt 



(-Ka + ie Q )8 ai + n5b* 2 + e^ 1 ^ k{o - 1) Ua^A 

(-K<j + ie )5a 2 + n8b\ - ie i4>1 y/ k(ct - 1) (Sa^ 

(-Ka + ie )Sb 1 + n6a* 2 + ie'^ 1 y/ n(cr - 1) (Sa^ 

(—kit + ieo)5b 2 + n8a\ — e -4 ^ 1 \J n(a — 1) ^6a^^j 



(in) 
(in) 



+ e Q Sa 2 - e 2 ^ 1 K(a - 1) 



2ft<5a™ 



(in) 



e 0< 5ai + e 2 ^ 1 k(<t - l)5h - V2^5a 2 n 
e Q Sb 2 + e~ 2 * k((j - l)5a 2 - V^Sbf 1 
eaSh - e- 2 ^ 1 n(a - l)S ai - V^5b 2 n 



r 



m 



where r stands for the cavity round-trip time, da] 71 cor- The evolution equations can be then expressed in matrix 
respond to the vacuum fluctuations entering the cavity form : 
through the coupling mirror. 

B. Quadrature fluctuations 



We define the amplitude and phase quadratures of a 
mode a with mean phase 4> a by respectively 



T ~dt 5J r = M5J + - iJWo" + 



^5J in (26) 



p a = Sue-**" + 5a* e 1 ^ 

q a = -i (Sae- 1 ^- -6a*e^') . 



(25) where 










eo 








-k(g - 2) 
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— eo 
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— KG 


-eo 


-n(e - 2) 











— eo 





eo 







-KIT 














-k(g - 2) 





— KG 


— eo 





eo 













eo 


— KG 


Co 





-k(<t - 2) 














CO 


— KG 


-co 













— eo 





Co 


— KG f 



M = 

\ 

6 J is the column vector of the quadrature components 
5 J = (Poi , q ai , Pa 2 , g Q2 , Pb! , 5b! , Pb 2 ,qb 2 ), (28) 
(5j7 m for the input pump fluctuations, 



(27) 



r 



and 8J ln the for the vacuum fluctuations entering the 
cavity through the coupling mirror 



5 7 in = (V™ a m v in a in 

n in n in n in n in 
Pb! ) Qbi 'Pb 2 J 1b 2 



(30) 



P { f\ q { f\p [ r\ti yn )- 



These differential equations are readily transformed 
(29) into algebraic equations by taking the Fourier transform. 
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In the frequency domain, the equations become 



2insj(n) = M'sj{ny 



where tl = 



n,. 



'-5JS n (ty+\l -SJ m (Cl). 

(31) 

is the noise frequency normalized to 
the cavity bandwidth il c = — and c = a is the normal- 
ized coupling constant. <jjf(f2), <5Jo m (fi) and SJ in (Q) 
are the Fourier transforms of the column vectors [2B1 EH1 
1201 The matrix M' is defined by 



differences. Thus, one can write two sets of equations: 
2iQSw±{Q) = M ± Sw ± (Q) (35) 



2- — -dw l ^ ± (n) + J-6wf(n) 



where Sw± are the column vector 

Sw-{n) = (r a (n),f fi (n),s a (a),s p (n)), (36) 



— C7 


— c 





c 








-cr-2 
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— £7 


— c 














— CT 
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—a 


— c 


-cr-2 











— c 
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— cr 





— a 














-a - 2 





— (J 


— c 





c 











— cr 
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— cr 


— c 





-a- 2 
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— a 


— c 





— <r 








— c 





c 


— cr 



(32) M+ = 



C. Correlations and anticorrelations 

In order to use the symmetry of the equations, one can 
introduce the symmetric and antisymmetric components 



Pa = -j= (Poi(O) +Pb 2 (ty) 



(33) 



and 



pp = -^(pa 2 (n)+p bl (o.)) 



(34) 



In a standard OPO, one expects intensity correlations 
and phase anti-correlations. Using our notations, this 
corresponds to having r a and as well as q a and qp 
squeezed [28|. This is also the case in a phase- locked 
OPO in the limit of small coupling p3 |. 

As previously, the equations verified by these quanti- 
ties can be expressed in matrix form. However, it is in- 
teresting to note that the sums are independent from the 



( -2(cr-l) 

-2(ct -1) 

c — c 

— c c 

( -2 -c -c\ 

-2 -c -c 

c c 

\ c c ) 



—c c 

c — c 

-2a 

-2ct / 



(37) 



6u)q 1 ± are the column vectors for the input pump fluctu- 
ations, 



5w%_ 



(0,0,0,0) 



(38) 



and 5w m the column vector for the input fluctuations 
entering the cavity through the coupling mirror 

5ts£(n) = (i£(n),j^(n),^(fi),^(fi)) 

cT<Dl"(fi) = (f^(Q),rf(n),s^(n),^(Q)) (39) 

These equations can be easily solved and one obtains 
the equations for the intracavity fluctuations. However, 
the relevant fluctuations are those outside the cavity. 
They can be easily calculated using the boundary con- 
dition on the output mirror: 

font =tf- rfi n w v 7 ^/ - fin (40) 
and the corresponding spectra are given by 

S/(fi) - (fout(W: ut (-n)) ■ (4i) 

The expressions are identical for the two sets of spectra 
^Pa,Qa^a-,s a and Sp q ?r „ jS ^ . 
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u = a, p. 



s g «(n) 
s«.(n) 



1 + 
i - 
i - 
i + 



i 

2(ft 2 + - l) 2 ' + 2 ((c 2 + o(o - l)) 2 + W- 2(c 2 - o-(c7 - l))ft 2 + ft 4 ) ' 

i (o- - 1) 2 + n 2 - c 2 

n 2 + a 2 ~ 2 ((c 2 + a {a - l)) 2 + ft 2 - 2(c 2 - a(a - l))ft 2 + ft 4 ) ' 



1 



ft 2 - 



2(1 + ft 2 ) 2(ft 2 + (ft 2 -c 2 ) 2 )' 

i i + ft 2 - c 2 



20 2 2 (ft 2 + (ft 2 - c 2 ) 2 ) ' 



(42) 



The behavior of the squeezed spectra, § r (<mt) and S^tcmt) , 
is shown on figure [21 For small coupling and low anal- 
ysis frequency, one observes that the spectra go to zero 
which indicates the existence of very large correlations. 
As usual, the spectra go to one as ft goes to infinity: 
the correlations exist only within the cavity bandwidth 
or conversely, only when the integration time is larger 
than the field life time in the cavity. The increase of c 
also results in a degradation of the correlations and anti- 
correlations: the wave-plate tends to mix the fields, and 
this phenomenon is increased as the coupling is increased. 
This phenomenon can be understood in terms of phase 
diffusion: as the coupling is increased, the phase diffu- 
sion is reduced which in turn reduces the intensity corre- 
lations. Finally, the phase anticorrelations (Fig. [2 left) 
are also degraded as the input pump power is increased 
: as ex is increased, the input pump noise becomes more 
and more important. 

This analysis has been performed by studying the fields 
generated by each crystal. Since several frequencies are 
involved, such a measurement would require the use of 
multiple local oscillators or of dephasing cavities [^j . It 
is thus easier experimentally as well as more fruitful for 
quantum information purposes to study the beams at 
each frequency, namely the couples (01,02) and (01,62), 
which are circularly polarized. 



IV. POLARIZATION ENTANGLEMENT 
GENERATION 

The non-classical properties of the two-crystal OPO 
are now analyzed in terms of polarization fluctuations. 
Polarization entanglement generation is demonstrated 
using various criteria introduced in the literature. 



A. Quantum Stokes parameters 

Before describing the properties of the system, let us 
briefly recall the quantum picture for polarization prop- 
erties. The polarization of light beams can be described 



in the general case by the Stokes parameters |3C 

"5*0 = Ix + ly 

S\ = I x — ly 

5*2 = ^+45° — ^—45° 

5*3 = I a + - la- ■ 



(43) 



where I x , y .+4:5 o .-45° ,0— describe the intensity along 
the x, y, +45°, —45°, er + , and o~ polarization. These 
four parameters arc not independent but are instead 
linked by the relation 



s$ = sf + s 2 2 + s 2 3 



(44) 



for a totally polarized field. 

The quantum polarization properties are defined via 
the quantum counter-parts of the Stokes parameters Q, 



2] 



So 


= ala x + 


&la y 






Si 


= ala x - 


al&y 






s 2 


= al& y + 


a y a x = 0+4500+450 - 




450 a_ 450 


s 3 


= i {ala x 


-al&y) =al +a<T+ - 




-K- 



(45) 



(Si, S*2, S3) verify commutation relations identical to 
those verified by orbital momentum operators: 



[Si , Sj 



2iS k 



(46) 



where i,j,k = 1,2,3 are cyclically interchangeable. 
These relations result in three uncertainty relations, 



A 2 SiA 2 S 2 > 
A 2 S 3 A 2 Si > 



(S3) 
{S2) 



, A 2 S 2 A 2 S 3 > 



(Si 



(47) 



where the notation A 2 (X) correspond to the variance of 
the operator X. 

B. Entanglement criteria 

In our case, Eq. 1191 shows that the output beams are 
circularly polarized, thus (S"' b ) = = (S^). It results 



8 



Stout) = S (o«.)((T = 1) 



5 9 ( »«,(<r = l.l) 




3.0 0.0 




FIG. 2: Normalized noise spectrum of the amplitude quadrature difference and phase quadrature sum for a = 1 (left) and 
a = 1.1 (right). Continuous curves corresponds to c = 1, long-dashed to c = 0.2 and short-dashed to c = 0. 



that the inequalities 1471 which contain variances S^' b are 
automatically verified. The only non-trivial Heisenberg 
inequality is then given by 



A 2 S^A 2 S 2 °'° > 



(48) 



where the notation A 2 (X) correspond to the variance of 
the operator X: A 2 {X) = ^(x - = (SX 2 ). 

Various entanglement criteria have been defined. We 
will restrict ourselves to three of them jjHl |32, HE 01 ■ 

The first two criteria arc, in the general case, sufficient 
(and not always necessary) conditions for entanglement. 
A general criterion is given by the product of two linear 
combinations of conjugate variables variances |33j (de- 
noted " product criterion" in the following) . In the case of 

I 



the Stokes operators, this inseparability criterion states 
that when 



A 2 (SI ± S b ).A 2 (S^ T S b 2 ) < 2 + | (S; 



(49) 



the state is entangled. The relevant Stokes parameters 
fluctuations, SS^' b and JS'^' 6 , can be expressed from the 
quadrature operators. In the case of circularly polarized 
beams, one has (u = a, b) 



\J\ 



-(Sq^m-Sq^Hm (50) 



It follows from these expressions that: 



ss+ = ssuci) + ss b (n) = (*js<?*>(fi) - 6p { ° ut \n)) ~ (sp^m - WrW) = V2(r^ ut) - r ( ; ut) ) 
ss; = ss%(n) - ss b (n) = -(Sq^(n) + 6qi° 2 ut \n)) + (n) + s^ffl) = V2(- q ^ + q ( ; ut) ) 



(51) 



Thus the relevant spectra can be calculated from the pre- 
viously calculated linear combinations of the amplitude 



and phase quadrature components. They arc given by: 



J 



s 5 +(n) 

S s -(fi) 



= 1 - 



n 2 - 



n 2 + (n 2 - c 2 ) 2 

{n 2 - c 2 ) + (a- l) 2 

[c 2 + a (a - l)] 2 + n 2 - 2 [c 2 - a(a - 1)] ft 2 + ft 4 ' 

I 



(52) 



Let us remark that the two expressions are equal for a = 1 and c = 0, i.e. for the uncoupled system op- 



FIG. 3: Criterion spectra for c = (left) and c = 1 (right). Continuous curves correspond to the sum criterion (|§ s +(51) + 
§ s -(f2)), long-dashed to the EPR criterion (S S a| S i,.S S a| S t,) and short-dashed to the product criterion (| ^S S +.S S -^). 



erated at threshold. One also remarks that S g +(f2) is 
independent of the pump amplitude, a. This is not sur- 
prising since this quantity is the difference of the intensity 
correlations between the signal and idler modes of each 
crystal. These correlations originates in the paramet- 
ric down-conversion phenomenon, independently of the 
pump power: the photons in the modes a\ and 62 («2 
and b\ respectively) arc emitted by pairs thus leading to 
intensity correlations between the two modes. These two 
quantities can be measured using standard methods for 
the measurement of quantum Stokes operators, see e.g. 

Using these expressions, one gets a simple expression 
for the product criterion 1491 

S s +(fi)S s -(n)<2. (53) 

The second criterion ("sum criterion"), which is widely 
used in continuous-variable entanglement characteriza- 
tion, is given by the sum of the above quantities |3ll |32| . 
It is in fact a special case of the previous criterion [33|. 
The criterion states that if 




then the two states are entangled. 

These two criteria are entanglement witnesses in the 
sense that they allow to specify whether the state is in- 
deed entangled. 

Finally, the last criterion is related to another partic- 
ularity of the system which may allow to make a joint 
QND measurement of the two non-commutating observ- 
ables Si and §2- It is shown in [34| that this is related to 
EPR- like "paradox". This criterion relies on conditional 
probabilities and can be written 

§Sf|Sf ■&s$\s!> ^ 1) 

where 



When this condition is verified, the correlations are suf- 
ficiently large that the information extracted from the 
measurement of the two Stokes operators of one field pro- 
vides values for the Stokes operators of the other which 
violate the Heisenberg inequality [35| . 

Figure 13 shows the various criteria as a function of the 
analysis frequency, f2 for a fixed pump power (a = 1) 
and for different values of the coupling. Very strong en- 
tanglement can be found for small couplings, c and small 
analysis frequencies, Q. As the coupling is increased, the 
entanglement is shifted to higher values of the analysis 
frequency. Let us note however that c = 1 corresponds 
usually to large angles (around 1° for typical output cou- 
plers): in the case of single crystal phase locked OPOs, 
stable phase-locking has been achieved with very small 
angles corresponding to c ~ 0. The entanglement is 
present for a wide range of the parameters thus show- 
ing the efficiency of the system. 



V. CONCLUSION 

We have presented and studied theoretically an orig- 
inal device based on an optical cavity containing two 
type-II phase-matched crystals as well as bircfrin- 
gent plates which add a linear coupling between the sig- 
nal and idler fields emitted in each crystal. The cou- 
pling induces a phase-locking respectively between the 
two signal fields and between the two idler fields. In 
this configuration, the system is found to directly gener- 
ate two-color polarization-entangled beams without the 
need for mode mixing. The two-crystal OPO would be 
thus a useful resource for the field of continuous- variable 
quantum information. Non-classical polarization states 
are well coupled to atomic ensembles which can be used 
as quantum memories and form the basic block required 
for quantum networks. Furthermore, we believe that the 
experimental implementation of this device is very inter- 
esting for potential parallel processing as it would be a 
new step towards extending the dimension of the phase 
space experimentally accessible. 
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